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Abstract. The Gaussian Effective Potential in a fixed transverse unitarity gauge is studied for the static 
three-dimensional U(l) scalar electrodynamics (Ginzburg-Landau phenomenological theory of supercon- 
ductivity) . In the broken-symmetry phase the mass of the electromagnetic field (inverse penetration depth) 
and the mass of the scalar field (inverse correlation length) are both determined by solution of the coupled 
variational equations. At variance with previous calculations, the choice of a fixed unitarity gauge prevents 
from the occurrence of any unphysical degree of freedom. The theory provides a nice interpolation of the 
experimental data when approaching the critical region, where the standard mean-field method is doomed 
to failure. 
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1 Introduction 

Since the discovery of high T c superconductors several 
unconventional models have been proposed in order to 
describe the unusual properties of cuprates. In particu- 
lar, the strong electron-electron coupling characterizing 
such materials requires new theoretical methods beyond 
the standard mean field approach. On the other hand 
the general phenomenology is still well described by the 
standard Anderson- Higgs mechanism: the supercurrent is 
carried by pairs of charged fermions whose non- vanishing 
expectation value breaks the gauge symmetry, thus en- 
dowing the gauge bosons with a mass. Thus the standard 
Ginzburg-Landau (GL) effective lagrangian still provides 
the best framework for a general description of the high-T c 
cuprates phenomenology. Moreover, as the GL action can 
be seen as a power expansion of the exact action around 
the critical point, the GL action must be recovered by any 
microscopic theory at least around the transition. Thus, 
regardless of the nature of the pairing mechanism, the GL 
action is a sound starting point for a general description 
of the high-T c materials. Of course we cannot trust the 
mean-field approach to the GL effective theory, and we ex- 
pect that many unconventional properties are connected 
with the breaking down of the simple mean field picture. 
Quite recently, the GL model has been extensively stud- 
ied both theoretically [T||2j and numerically [Hj in order to 
clarify the universality class of the superconducting transi- 
tion and the role of the critical fluctuations in the high-T c 
cuprates, as well as the order of the transition itself 0]. 

Actually the high-T c cuprate superconductors are char- 
acterized by a very small correlation length £ which al- 



lows the experimentalists to get closer to the critical point 
where the thermal fluctuations cannot be neglected and 
the mean field approximation is doomed to fail As far 
as we know there is no full evidence that the universal 
critical behaviour has been reached in any real sample |S] , 
but it is out of doubt that an intermediate range of tem- 
perature is now accessible, where thermal fluctuations are 
not negligible even if the sample is still out of the truly 
critical regime. Thus, in order to describe some unconven- 
tional properties of the high-T c superconductors, we need 
to incorporate the role of thermal fluctuations, but unfor- 
tunately we cannot rely on the standard renormalization 
group methods which would only describe the universal 
limiting behaviour that could not be observed yet in any 
real sample. We need some kind of interpolation scheme 
for the non-universal regime where the behaviour depends 
on the physical parameters of the sample, and we would 
prefer a non-perturbative method in order to deal with 
any strong coupling. 

In this paper we study the Gaussian fluctuations by 
means of a variational method, the Gaussian Effective 
Potential (GEP), which has been discussed by several au- 
thors as a tool for describing the breaking of symmetry in a 
simple scalar theory [7||H] • As a toy model for electro- weak 
interactions, the scalar electrodynamics in four dimensions 
has been studied by Ibanez-Meier et al. [§] who computed 
the GEP by use of a general covariant gauge. However 
their method gives rise to an unphysical, and undesirable, 
degree of freedom. 

We compute the GEP for the U(l) scalar electrody- 
namics in three space dimensions where it represents the 
standard static GL effective model of superconductivity. 
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In order to make evident the physical content of the the- 
ory, thus avoiding the presence of unphysical degrees of 
freedom, we work in unitarity gauge. This has been shown 
to be formally equivalent to a full gauge-invariant method 
once all the gauge degrees of freedom have been integrated 
out 10' . The variational method provides a way to eval- 
uate both the correlation length £ and the penetration 
depth t as a solution of coupled equations. The GL pa- 
rameter kgl = is found to be strongly temperature 
dependent in contrast to the simple mean-field descrip- 
tion On the other hand the model predictions are 
in perfect agreement with some recent experimental data 
[H] . which can be nicely interpolated by our variational 
calculation. 

The comparison with the experimental data is of spe- 
cial importance as it provides a test for the GEP vari- 
ational method itself. The predictions of the method in 
3+1 dimensions, in the context of electro- weak interac- 
tions, have been discussed by several authors |7|5pl3pi4| . 
but no real comparison with experimental data will be 
achievable until the detection of the Higgs boson. Thus 
high-T c cuprate superconductors represent our next best 
choice in order to test the reliability of the method. 

The paper is organized as follows: in section the GL 
action and partition function are introduced and discussed 
within the unitarity gauge; in section|3the Gaussian vari- 
ational method is applied to the three-dimensional GL 
model, and the coupled variational equations are derived; 
finally, in section 01 the phenomenological predictions of 
the method are compared with some experimental data 
for the GL parameter. 



define two real fields p and 7 according to <fi — pe n . The 
unitarity gauge is recovered by the gauge transformation 



A A - — Vy(x) 



(5) 



and the original transverse vector field Aj_ acquires a lon- 
gitudinal component proportional to V7. Thus the 
original measure in Eq. becomes 

J D[4>,<f>*,A ± ] - J D[j}pD[p}D[A ± } 

-» const x J pL>[ / 9]L»[A L ]L>[A ± ] (6) 
In unitarity gauge the action, Eq. l(T|l. now reads 
S = J d 3 x j^(Vp) 2 + \m\p 2 + \ B p" 

+i e |p 2 (A| + Ai) + i(VxA ± ) 2 | (7) 

and the longitudinal field A L may be integrated out ex- 
actly yielding a constant factor and an extra 1/p factor 
for the measure Finally, dropping the constant fac- 
tors, the partition function may be written as 



Z = I D[p, AJcxp 
1 



I( V p) 2 + im|p 2 



+ XBP 4 + ^e 2 B p 2 A 2 ± + -(\7 x A ± f 



(8) 



2 The GL action in unitarity gauge 

Let us consider the standard static GL action J5| 



S = 



\F^Fi* + ±(D l >4>r(D»4>) 



+ -m|0*0 + A B (0*0) 2 



(1) 



Here is a complex (charged) scalar field, whose covariant 
derivative is defined as 



D„ 



dn + ie B A, 



(2) 



and /i, v — 1,2,3 run over the three space dimensions. The 
magnetic field components F av — d^A v — d v A a satisfy 



\FuuF^ = |V x A| 



(3) 



and the partition function is defined by the functional in 
tegral 



D< 



(4) 



We may assume a transverse gauge V • A = 0, and then 
switch to unitarity gauge in order to make <p real. Let us 



We may enforce the transversal condition on the vector 
field by a gauge fixing term in the action and, restoring 
p = (f), the action reads 



S = 



d 3 x 



1 2 
2 % 



+ A B 4 + l-ej 



'■A 1 



i(VxA) 2 + l(V-A) 2 



(9) 



The partition function is now expressed as a functional 
integral over the real scalar field and the generic three- 
dimensional vector field A, with the extra prescription 
that the parameter e is set to zero at the end of the cal- 
culation. Inserting a source term we may write 

Z\j] = J D[cj>,A a }expl.-S + J d 3 xj<f^ (10) 

with S given by Eq. 10. The free energy (effective poten- 
tial) follows by the Legendre transformation 



F[tp\ =-lnZ + / d 3 xjip 



(11) 



where ip is the average value of <$> in presence of the source 
j. The superconducting phase is characterized by an ab- 
solute minimum of T for ip =/= 0. 
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3 The Gaussian method 

From the partition function l|10[l . the GEP may be evalu- 
ated by the S expansion method discussed in Refs. [§1TH$] . 
In our case the GEP represents a variational estimate of 
the free energy l(TT|) . 

First we introduce a shifted field 

= 0-^ (12) 

then we split the lagrangian into two parts 

C = Co + Hint (13) 

where Co is the sum of two free-field terms describing a 
vector fi< 
mass fl: 



vector field with mass A and a real scalar field </> with 



Co = 



-~(Vx A) 2 + l -A 2 A^ + Uv-A) 2 



2c 



1 ~ 7 ' 2 ■ -fi 2 s 2 



5 CW 



The interaction then reads 



(14) 



Ant = «0 + Vx<f> + V2(f> +v 3 cb 
+± (e%tp 2 - ^) 
1 



where 



m 2 B tp 2 + X B tp A 



vi = m 2 B ip + i\ B tp 3 

v ? = -m| + 6A B <y9 2 - ^fl 2 

«3 = 4A B <p 
Vi — As- 



(15) 

(16a) 

(16b) 

(16c) 

(16d) 
(16e) 



We now expand the free energy to first order in £j nt 
following standard perturbation theory procedures. We 
obtain 

T[tp] = \ Trln [g-\x, y)] + ~ Trln [G^(x, y)] 
+ / d 3 x {v + v 2 g(x, x) + 3v i g 2 (x, x) 



+ -e 2 B (g(x 1 x) + p 2 -A 2 )G^(x,x)j (17) 

where g(x,y) is the free-particle propagator for the scalar 
field, and G tlv (x, y) is the frec-partiele propagator for the 
vector field 



Gj,{x,y) = 



+ ( - - 1 ) k^k v 



(18) 



In the limit e — > 0, up to an additive constant 
Tr In [G;J (s, y)] = 2V y — ln(fc 2 + ^ 



) (19) 



where V is the total volume. Dropping all constant terms, 
the free energy density V c s — J-/V (effective potential) 
reads 



V eS [ V ]=h{Q) + 2h{A) 
1 



\ B p 



1 2 

m B tp 



fi 2 



+ nx B p 2 + &\ B io{n)} h{n)] 

+ (e 2 B tp 2 + e%I (f2) - A 2 ) I Q (A) (20) 
where the divergent integrals /„ arc defined according to 

d 3 k 1 



h(M) = 



(2tt) 3 M 2 + k 2 



h{M) = \J ^b,(M 2 



(27T) S 



k 2 ) 



(21) 



(22) 



and are regularized by insertion of a finite cut-off A. 

The free energy (1201) now depends on the mass param- 
eters fl and A. Since none of them was present in the 
original GL action of Eq. © , the free energy should not 
depend on them, and the minimum sensitivity method |16| 
can be adopted in order to fix the masses: the free energy 
is required to be stationary for variations of fl and A. On 
the other hand the stationary point can be shown to be 
a minimum for the free energy and the method is equiv- 
alent to a pure variational method UJ. At the stationary 
point the masses give the inverse correlation lengths for 
the fields, the so called coherence length £ = X/Q and 
penetration depth I — 1/A. 

The stationary conditions 



on 2 

dV eS _ 
OA 2 

give two coupled gap equations: 







(23a) 
(23b) 



fl 2 = 12X B I (O) +m 2 B + l2X B p 2 + 2e 2 B I (A) (24a) 



a1 2 2 

A = e B tp 



,J (ft) 



(24b) 



For any tp value, Eqs. (|24(l must be solved numerically, 
and the minimum point values fl and A must be inserted 
back into Eq. H2()(l in order to get the gaussian free en- 
ergy V e ff (tp) as a function of the order parameter tp. For 
a negative and small enough m B , we find that V e e has a 
minimum at a non zero value of tp = p lri \ n > 0, thus indi- 
cating that the system is in the broken-symmetry super- 
conducting phase. Of course the masses fl, A only take 
their physical value at the minimum of the free energy 
Pmin- This point may be found by requiring that 



dtp 2 







(25) 
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Fig. 1. The GL parameter according to Eq. 127^ for kq — 100, 
Co = 1-36 nm, T c = 121.5 K and A£o = 57 (full line). The circles 
are the experimental data of Ref. for Tl2Ca2Ba2Cu30io- 



where as usual the partial derivative is allowed as far as 
the gap equations l|24|l are satisfied 'T. The condition 125( 1 
combined with the gap equation 124a(l yields the very sim- 
ple result 

Q 2 

(26) 



Vmin 



8A, 



However we notice that here the mass Q must be found 
by solution of the coupled gap equations. Thus Eqs. i|26(l 
and (|24J) must be regarded as a set of coupled equations 
and must be solved together in order to find the physical 
values for the correlation lengths and the order parameter. 

Insertion of Eq. (|26J1 into Eq. (|24b(l yields a simple 
relation for the GL parameter kql 



"GL 



1 



w 1 . 

r mm 



(27) 



where kq = eg/(8As) is the mean-field GL parameter 
which does not depend on temperature. Equation l|27|l 
shows that the GL parameter is predicted to be temper- 
ature dependent through the non trivial dependence of fi 
and ^min- At low temperature, where the order parameter 
Vmin is large, the deviation from the mean-field value Ko 
is negligible. Conversely, close to the critical point, where 
the order parameter is vanishing, the correction factor in 
Eq. (|27|l becomes very important. 



4 Comparison with experimental data 

In order to make contact with the phenomenology of the 
high-T c cuprate superconductors we need to fix the bare 
parameters of the GL action. The standard derivation JS] 
of the GL action (JTJ from a microscopic model gives a di- 
rect connection between microscopic first-principle quanti- 
ties and phenomenological bare parameters. The bare cou- 
pling e b turns out to be related to the elementary charge 



of fermions, to the critical temperature T c and to the zero 
temperature coherence length £ through |15j 



^y/k B T c £ . 
nc 



(28) 



The other parameters may be fixed by knowledge of some 
zero temperature phenomenological quantities like coher- 
ence length and penetration depth. The bare mass param- 
eter m B may be regarded as a linear function of temper- 
ature 

m 2 B =m 2 c +(l- !j (m§ - m 2 c ) (29) 

where m,Q is the value which is required in order to find 
Q = l/£o from the gap equation l|24a(l at (p = ip m i n , and 
m 2 is the value of m 2 B at the transition point. In mean-field 
approximation m 2 = 0, while here the fluctuations shift 
the transition point to a non vanishing m 2 B value. The 
bare couplin g Ag is regarded as a constant, and is fixed 
through Eq. (|24b(l by requiring that at zero temperature 
(i.e., for m 2 B = to§) the penetration depth Iq = 1/A. This 
way the method provides a one-parameter interpolation 
scheme for the superconducting properties, as the cut-off 
A still has to be fixed. The cut-off A is a characteristic 
energy scale of the sample, and may be determined by a 
direct fit of the experimental data. 

Unfortunately there are not too many available exper- 
imental data on the behaviour of superconductors close 
to the critical point. Even for the high-T c cuprate su- 
perconductors the measurement of both coherence length 
and penetration depth up to the pre-critical region is not 
an easy task. The GL parameter has been reported by 
Brandstatter et al. |12j as a function of temperature for 
the high T c material Tl 2 Ca 2 Ba2Cu30io (T c = 121.5 K, 
kq = 100, £o = 1-36 nm). These data are shown in Fig.^ 
together with the interpolation curve obtained by Eq. (|27|l 
for A£o = 57. We observe that the experimental GL pa- 
rameter is almost constant kql ~ up to T — 0.8T C 
where it starts decreasing. Thus for T/T c > 0.8 the mean- 
field approximation breaks down and the thermal fluctu- 
ations become important. On the other hand Eq. I|27(l fits 
the experimental data up to 98% of the critical temper- 
ature. Beyond that point there are neither experimental 
data nor reliable theoretical predictions. As we approach 
the critical point some universal behaviour should be ex- 
pected and the role of thermal fluctuations becomes too 
important to be dealt with by the present method |17j . 

In perturbative calculations, at one-loop order, the uni- 
tarity gauge is known |18j to give rise to wrong predictions 
around the critical point. Thus it could be argued that, 
because of the bad ultraviolet behaviour of the gauge prop- 
agator, higher order diagrams should be included to cancel 
the divergences even in the present one loop calculation. 
However, at variance with perturbative approximations, 
variational calculations are known to be less sensitive to 
higher order corrections. The GEP provides very sensi- 
ble results even when the one loop effective potential fails 
entirely [7j. Moreover, while the way of dealing with di- 
vergences is crucial at the critical point, the present cal- 
culation deals with a pre-critical region where the cut-off 
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regulator is finite and plays the role of a physical length 
scale which is naturally determined by the structure of 
the condensed matter system. Thus divergences are not 
an issue: our unitarity gauge lagrangian must be regarded 
as regularized and is known flB] to give the same results 
predicted by other gauge choices. 

It must be kept in mind that all variational methods 
are quite sensitive to the choice of the trial functional. The 
Gaussian functional is not by itself gauge invariant, and a 
gauge dependent result is expected anyway by this varia- 
tional method. The question of determining the best gauge 
choice has been addressed by Ibahez-Meier et al. [H] who 
found the Landau gauge (d ■ A) = to be optimal in four 
dimensions by variational arguments. In three dimensions 
that is equivalent to the transverse gauge we started with 
in section [21 Unfortunately the Gaussian functional also 
breaks the U(l) symmetry of the lagrangian, and unphys- 
ical massive degrees of freedom show up unless the uni- 
tarity gauge is chosen. Thus our gauge prescription must 
be regarded as a compromise which allows to make con- 
tact with the real phenomenology. The comparison with 
experimental data is encouraging in this respect. 

As the GEP provides a nice way to interpolate the ex- 
perimental data beyond the mean-field regime, we expect 
the method to be reliable for the description of symmetry 
breaking in 3 + 1 dimensions where the scalar electrody- 
namics may be regarded as a toy model for the standard 
electro-weak theory. Since there is no reason to believe 
that the real universe is very close to the transition point, 
the Gaussian method may be regarded as a valid tool for 
describing the effects of fluctuations beyond mean-field 
approximation. 
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